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1. Introduction 

Galaxies are not distributed randomly but are concentrated within elon- 
gated filamentary chains which consist of groups and clusters of galaxies; 
the space between filaments is devoid of galaxies. Such distribution can be 
called cellular (Joeveer and Einasto, 1978; Zeldovich et al. , 1982): a cell is 
a large low-density region surrounded by superclusters. Examples are the 
Northern Local Void surrounded by the Local, Coma and Hercules super- 
clusters (Lindner et al. , 1995), and the Bootes void (Kirshner et al. , 1981) 
surrounded by the Hercules and Bootes superclusters. Superclusters and 
voids form a continuous network of alternating high- and low-density re- 
gions; the mean diameter of voids between clusters in the supercluster-void 
network is about 100 h^^Q Mpc (Zeldovich et al. , 1982). 

It is not clear whether superclusters and voids form a regular or irregu- 
lar network. According to the classical paradigm of the formation of large 
scale structure the distribution of density waves is Gaussian on all scales, 
and thus the supercluster-void network should have a random character. 
The observed network was formed by density waves of a wavelength range 
corresponding to the scale of the network. Therefore, it was a great surprise 
when it was found that the distribution of high-density regions along pen- 
cil beams around the northern and southern Galactic pole is fairly regular: 
high- and low-density alternate with rather constant step of 128 Ii^qq Mpc 
(Broadhurst et al. , 1990). The regularity of the structure is so far well 
established only in the direction of Galactic polar caps, while in other di- 
rections the regularity is much less pronounced. In order to find the de- 
gree of the global regularity of the supercluster-void network 3-dimensional 
data of the distribution of high-density regions are needed. For this pur- 
pose Abell-ACO clusters of galaxies (Abell, 1958; Abell et al. , 1989) can 
be used. These clusters cover the whole celestial sphere outside the Milky 
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Figure 1. Distribution of clusters in high-density regions in supergalactic coordinates 
(Einasto 1995). Only clusters in superclusters with at least 8 members are plotted. The 
supergalactic Y = Q plane coincides almost exactly with the Galactic equatorial plane; 
the Galactic zone of avoidance is marked by dashed lines. 

Way zone of avoidance. Here I give a summary of our principal results, a 
more detailed analysis is published elsewhere. 

2. Distribution of clusters of galaxies 

The distribution of clusters of galaxies located in very rich superclusters 
with at least 8 member-clusters is shown in Figure 1 (Einasto et al. , 1994; 
Einasto, 1995; Einasto et al. , 1997b), see also (Tully et al. , 1992). This 
Figure shows clearly that high-density regions are separated from each other 
by a fairly constant intervals of ^ 120 K^qq Mpc; in other words, they form 
a quasi-regular network of superclusters and voids. 

We supplemented this qualitative description by a quantitative analysis, 
using the power spectrum of clusters of galaxies (Einasto et al. , 1997a). 
On short wavelengths the spectrum can be approximated by a power law 
with index n = —1.8. On wavenumbcr = 0.05 /iioo Mpc^"*^ it has a 
sharp peak; this wavenumber corresponds to the wavelength Aq = 27r/fco = 
120 K^QQ Mpc. On longer scales the error corridor of the spectrum is large 
and here the spectrum is compatible with the Harrison-Zeldovich spectrum 
of power index n = 1. 

The presence of a sharp maximum in the power spectrum of matter is 
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the main finding of our study of the distribution of clusters of galaxies. Our 
result has found independent support from other data (Landy et al. , 1996; 
Gaztanaga and Baugh, 1997; Peacock, 1997; Rctzlaff, 1997). Comparison 
with simple toy models shows that a peaked power spectrum is possible only 
if high-density regions form a quasiregular rectangular network (Einasto et 
al. , 1997d; Einasto et al. , 1997c). In this case the correlation function 
of objects located in high-density regions is oscillating. Evidence for an 
oscillating cluster correlation function has been accumulating already for 
some time (Kopylov et al. , 1988; Mo et al. , 1992; Einasto and Gramann, 
1993; Fetisova et al. , 1993). 

3. Comparison with CMB data 

The angular power spectrum of the cosmic microwave background (CMB) 
has been measured by a number of teams. We compared the CMB spectrum 
with optical data using three models for the power spectrum: (a) a scale- 
free initial spectrum, (b) a double power law approximation to the cluster 
spectrum, and (c) a spectrum based on the observed cluster spectrum. For 
a set of cosmological parameters, we calculate the matter transfer function, 
and the matter and radiation power spectra for all three models (Atrio- 
Barandela et al. , 1997). We assume the Universe has a flat geometry. In 
calculations we used the package CMBFAST (Seljak and Zaldarriaga, 1996). 
To estimate the goodness of a particular set of cosmological parameters we 
calculate the parameter for all three principal models, see Figure 2. 

In Figure 3 we compare matter power spectra and temperature anisotropy 
spectra for our three basic models with the data. The cosmological parame- 
ters were chosen to reproduce the CMB data (Netterfield et al. , 1997). We 
see that the temperature anisotropy spectra are very similar in the range 
of multipolcs observed in Saskatoon. In other words, the present CMB data 
are not sufficient to discriminate between models. The matter power spectra 
are also similar on short wavelengths but on medium and long scales they 
are very different. The scale-free model with large cosmological constant has 
a broad maximum at large wavenumber (k ~ 0.01 /iioo Mpc~^); the maxi- 
mum of the first acoustic oscillation occurs at A; « 0.1 /iioo Mpc~^ and is of 
rather small relative amplitude. Both scales are outside the allowed range of 
the observed spike in the cluster spectrum: ko = 0.052 it 0.005 hioo Mpc~-^ 
(Einasto et al. , 1997a). No combination of cosmological parameters can 
reproduce the spike at k = ko: the existence of a broad maximum is an 
intrinsic property of all scale-free models. Thus the observed spike is not 
related to acoustic oscillations in the baryon-photon plasma as assumed 
by Szalay (1997) but must have a different origin. On the other hand, the 
cluster and double power law spectra fit the observed cluster spectrum by 
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Figure 2. Goodness-of-fit contours of at 68% and 95% confidence level. The 
statistics measures the deviation of the expected temperature anisotropy amplitude of 
a given model from the Saskatoon data. The first row displays the results for the scale 
free model; the second row for the double power law model; and the lowest row for the 
cluster spectrum based model. In the first column we plot models with varying Ifubble 
constant and baryon fraction for a spectral index n = 1 at large scales and no cosmological 
constant. In the middle column the same diagrams were repeated for n = 1.2. Dashed 
lines indicate the nucleosynthesis bounds. The last column displays the results for models 
with different values of the cosmological constant. On all these models, the age of the 
Universe was chosen to be 14 Gyr. 



construction and reproduce the CMB data, i.e they fit equally well both 
datasets. 

Thus the comparison of optical and CMB data brings us to the conclu- 
sion that CMB data are not in conflict with the presence of a spike in the 
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Figure 3. Comparison of matter power spectra and radiation temperature anisotropies 
with cluster and CMB data. Dots with la error bars give the observations; the measured 
cluster spectrum (Einasto et al. 1997a) in the left panel and the Saskatoon data on CMB 
temperature anisotropies (Netterfield et al. 1997) in the right panel. The scale-free model 
spectra (short-dashed lines) were computed using the following parameters: h = 0.6, 
Q.b = 0.07, Q.C = 0.23, and 0,^ = 0.7. The cluster spectrum (solid lines) was calculated 
using h = 0.6, Slf, = 0.08, Oc — 0.92, and Oa = 0; and the double power law models 
(long-dashed) using h — 0.6, Qb = 0.05, Sic = 0.95, and Q,a ~ 0. To compare matter 
power spectra and observations we used a bias factor bd « 3 for cluster spectrum. 

matter power spectrum, and that the present combined cluster and CMB 
data favour models with a built-in scale in the initial spectrum. We repeat 
that a regular supercluster-void network can be formed only by a power 
spectrum with a sharp maximum (Einasto et al. , 1997d). As noted by Sza- 
lay (1997) correlated phases are also crucial to form a regular network of 
high- and low-density regions. 

Double inflation models provide a possible scenario where the formation 
of a spike could have taken place. One version of a double inflation model is 
suggested by Starobinsky (1992). This model produces a spectrum rather 
similar to the initial spectrum found from data (Atrio-Barandela et al. , 
1997). The study of the distribution of matter on large scales is of crucial 
importance since it could provide a direct test of more complicated models 
of inflation. 

I thank H. Andernach, F. Atrio-Barandela, M. Einasto, S. Gottlober, V. 
Miiller, A. Starobinsky, E. Tago and D. Tucker for fruitful collaboration and 
permission to use our joint results in this review article, and A. Szalay for 
discussion. This study was supported by the Estonian Science Foundation. 
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